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NASA DISCLAIMER 


The conclusions stated herein are those of the Contractor and 
are not necessarily those of NASA. They are being published to 
direct attention to the problem of the effect of heavy rain on 
the aerodynamic performance of an aircraft. The theory proposed 
herein incorporates certain assumptions and extrapolations because 
suitable data do not exist. Because of this, the results and con- 
clusions reported herein are in question. They are published, 
however, in the hope that other researchers will be inspired to 
suggest and try new theoretical approaches or experimental programs 
to obtain the needed verifications. 
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FOREWORD 

The research reported in this document is concerned with 
the aerodynamic penalties of very heavy rain. Although wind 
shear has been implicated as the cause of a number of aircraft 
crashes in landing configuration , in most cases extremely heavy 
rain was also present. In the accident reconstruct ions / no 
accounting for the effect of heavy rain was made and for this 
reason derived wind shears may be too large. This research is 
concerned with the frequency and intensity of very heavy rains 
and their effect on landing aircraft. The effects consist 
largely of momentum losses by the aircraft to the rain/ increased 
drag due to roughening of the aircraft surface and decreased lift 
at high angles of attack due to a roughened surface. The magni> 
tude of these effects was established by analysis of raindrop 
trajectories in potential flow about an aircraft and the runback 
of those drops impinging on the aircraft. Momentum and drag 
penalties were introduced into a landing simulation model. Their 
effect was assessed in terms of the departure of the landing 
touchdown point from the touchdown point for the case of no rain. 
The departures are also compared to those previously calculated 
for landings done into severe low-level wind shear. 

The research was conducted by the University of Dayton 
Research Institute for the National Aeronautics and Space 
Administration/ Wallops Plight Center/ Wallops Island/ Virginia/ 
under the technical direction of Nr. W. E. Nelson and R. E. Carr 
of the Engineering Division/ Suborbital Projects and Operations 
Directorate. For the support of this effort/ the authors are 
indebted to Nr. A. Richard Tobiason of the Office of Aeronautics 
and Space Technology (OAST)/ NASA Headquarters/ Washington/ D.C. 
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SECTION 1 
EXORDIUM 

In recent years, wind shears associated with strong thunder- 
storm downdrafts have been implicated as the cause of several 
aircraft accidents. The Eastern Flight 066 accident at Kennedy 
Airport (NTSB, 1976) is a prime example. In the National 
Transportation Safety Board's reconstruction of the flight 
recorder data from the accident, extraordinarily large wind 
shears v^re estimated. The reconstruction considered no other 
external factors besides the wind. The performance degradation 
due to the heavy rain cell experienced by Eastern 066 was not 
taken into account. We feel it possible that the derived wind 
shears are coo large because the effect of the very heavy down- 
pour was ignored. 

An extensive literature search revealed only one other 
investigation (Rhode, 1941) which considered the effect of heavy 
rain on aircraft performance. That investigation dealt with the 
case of an aircraft encountering heavy rain at moderate altitude 
(about 5000'). It concluded that although heavy rain has a 
significant effect, its exposure time is insufficient to force 
the aircraft to the ground. An aircraft in landing con- 
figuration, however, does not have a wide margin of performance 
i:^ which to overcome the aerodynamic penalties due to heavy rain. 
Thus, in this report we consider the importance of heavy rain to 
landing aircraft. 

Rain can affect an aircraft in at least four ways: (a) 

raindrops striking the aircraft impart a downward and backward 
momentum; (b) a thin water film results from the rain increasing 
che aircraft mass; (c) the water film can be roughened by drop 
impacts and surface stresses producing aerodynamic lift and drag 
penalties; (d) depending on aircraft orientation, raindrops 
striking the aircraft unevenly impart a pitching movement. In 
this investigation, we used recent advances in computational 














fluid dyn.«lc8 to assess the first three penalties. In a 
following report the role of rain In several aircraft accidents 

will be addressed. 
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SECTION 2 

CHARACTERISTICS OP VERY HEAVY RAIN 

It is first necessary to establish the nature and frequency 
of very heavy rain. Only short duration downpours associated 
with convective cells possess the intensity to seriously affect 
aircraft performance. We have chosen to cateqorize such rains by 
their rainfall rates in millimeters per hour. 

Several authors have attempted to establish the frequency 
of occurrence of rainfall rates. Jones and Simms (1978) par- 
formed an analysis of one-and four-minute rainfall rates for many 
O.S. stations. The data used was only for a one-year period 
however, and the highest rate observed during that year was 238 
mm/hr at Miami, Florida. Hershfield (1972) performed an analysis 
to obtain the expected mean maximum 5-minute rainfall rate for 
stations throughout the U.S. He also estimated the maximum rain- 
fall rates over 1-minute periods to be about 50% higher than 
those over 5-minute periods, ■niese results give yearly mean 
maximum 1-minute rainfalls of from 150 to 250 mm/hr in the 
eastern United States. Information about rainfall rates for 
periods shorter than one minute is unavailable* It is con- 
ceivable that very short period (20 to 30 second) rates are even 
greater than the maximum one-minute rates. Thus, the world 
record rainfall rate of 1828.8 nm/hr or 73.8 in/hr at Unionville, 
Maryland, (Riorden, 1970) suggests an upper limit of 2000 m/hr. 

A 2000 m/hr rain will hereafter be cha:racterized as incre- 
dible. While a 500 m/hr rain rate is an unlilcely occurrence at 
a single station, the chances that it will occur at one or more 
stations in the eastern United States each year are considerably 
greater. Such a rate will be termed severe. A rainfall rate of 
100 m/hr or greater is expected at least once a year in most of 
the eastern United States. We term this rata as heavy. 
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SECTION 3 

SIZE DISTRIBUTION OP WATER DROPLETS 


To analyze the effect of heavy rain on aircraft perfor- 
mance, the size distribution of water droplets under different 
rainfall rates must first be established. The classical paper on 
the size distribution of rain drops is that of Marshall and 
Palmer (1948). Their results show that the size distribution can 
be approximated by the exponential function. 


dN(D) - N e"**^ 

Sd- ® 


( 1 ) 


where ♦ • 41 R number of drops within diameter 

range dD, N© is an empirical constant (.0), D is the drop 

diameter, and R is the rainfall rate in mm/hr. The drop distri- 
butions derived from Eq. (1) for several different rainfall rates 
are presented in Figure 1. The distributions described by (1) 
were originally derived from extratropical rains but Nercerec 
(1975) found it valid for tropical showers as well. 


The cerminal velocity of raindrops of drop diameter, D, has 
been established by Markowitz (1976) 

V(D) - ».5« U-.XP 

wh.r. V(0) ii th. terminal velocity. A correction Cor terminal 
velocity aloft is given in Markowitz 

0.4 

V(D) . V^(D) (^) 

where V^(D) is the terminal velocity for density p , p is the 
density aloft, and p^ is the density at the surface. Equation (3) 
allows terminal velocity adjustment for aircraft operating at 
higher flight levels. 


It is necessary to )cnow the percentage that each size 
droplet comprises of the total rain volume. The reason for this 


% 
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will b«con« evident later in the impingenent efficiency section. 
The percentage can be calculated given the terminal velocity of 
each size droplet and the number of droplets in each size range. 
Following Markowitz (1976) the frequency distribution for the 
fraction of total rain volume due to raindrops of diameter D is 

M(D)dD - N(D)|»(D/2)^ V(D)dD /*N( D)^*( D/2) ^ V(D)dD (4) 

where M(D) is the percent volume of total rain volume of drops of 
diameter D. Figure 2 shows the percentage volume contributions 
foi three rainfall rates. 


100 inn/hr 
500 mm/ hr 
1000 mm/hr 
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SECTION 4 

IMPINGEMENT EFFICIENCY 

Obviously, not all drops in the path o£ an aircraft strike 
it* Some, especially smaller drops, are carried over the 
aircraft by the flow of air caused by the aircraft motion. A 
first step in calculating the aircraft momentum loss due to rain 
is calculating the ratio of the rain that strikes the aircraft to 
the total that would strike the aircraft in the absence of the 
aircraft induced airflow. This ratio is called collection or 
impingement efficiency. 

We first calculated impingement efficiencies for a range of 
drop diameters from *5rom to 8rom. By knowing the relative percen- 
tage of each drop size, it was possible to establish impingement 
efficiencies by rainfall rate. This was done by summing the pro- 
duct of impingement efficiency and relative volume percentage for 
all drop sizes. 

A water drop trajectory program together with a potential 
flow model were used to calculate local and overall collection 
efficiencies of water drops by the airfoil of an aircraft. The 
potential flow model calculated airflow about the airfoil. The 
airflow was used in calculating the paths of raindrops. The 
calculations were done in proximity to the aircraft for a range 
of drop size diameters. For momentum calculations, a summation 
of the products of collection efficiency times percentage contri- 
bution by drop diameter was made for each rainfall rate. 

The Hess (1972) potential flow model was used because of 
the ease with which off-body velocities were obtained. This 
program has the ability to calculate the airflow about an 
arbitrary 3-D shape at up to 200 off-body points. These points 
are arranged in a two-dimensional grid so that airfoil sections 
are centered on the grid. The arrangement is shown in Figure 3* 



Figure 3. Two-Dimensional Grid of Off-Body Points for a Wing Cross-Section 











An arbitrary aircraft configuration is specified by geometry 
input which includes fuselage length and breadth, wing chord, and 
wing thicknesses. The angle of attack and speed may also be set 

w 

arbitrarily. 


Raindrops are assumed to be spherical except for very large 
drops for which an empirical correction for non-sphericity is 
applied. Their motion is governed by Newton's second law. Under 
the assumptions that forces acting on a drop arise solely from 
fluid drag and gravity and that the drop is small enough so as 
not to disturb the airflow field, the second law can be written 
as 


where 


n, 

”d —I 

dt 


^Drag “ “d ^ 


m. 


drop mass 


(5) 


^drag “ 

3^ « drop velocity 

g » gravity 

which expressed in non-dimensional form is (Bergrun, 1951; 
Norment, 1976) 

d3. C^Re 

- 5iiT ‘9. - 

Where Re * °^^a"^d^^free^“a Reynold's number, 

K « is the inertial parameter, 

F ■ ^/Lg is the Froude number, 

N free ^ 


t'V 


free 


free 


'free 


p^ the air density. 
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p^ the water density, 

U the air velocity at the drop position, 

A 

V. the free stream velocity, 
tree 

R. the drop radius, 

U 

L a characteristic length either the wing chord or fuselage 
length, 

and V is the viscosity. 

The drag coefficient, C^, is a function of the drop 
Reynolds's number (Walsh, 1977) 

Cp ■ (1 + 0.1 5Re ). 

For large Reynolds number, i.e., Re>200, the drag must be 
corrected for non-sphericity of the drops. is iterated in the 
empirical relation (Berry and Pranger, 1974) 

2 3 

InRe-ajQ + a2j, InXj + a22 (lnX2) + a23(lnX2) 
where X 2 « CjjRe^ a 2 Q » -0.236534 x 10^, a 2 j^ » 0.767787, 

a 22 - 0.535826 x lO"^, and a 23 - 0.763554 x lO"^ 

until a solution is obtained. 

Eq. (6) is integrated numerically by a 4th order 
Runge-Kutta scheme to obtain 0^(t’fat) and drop positions are 
incremented by (AcJcley and Templeton, 1979), 

Z. - Z. * M.dt 
d Q d 

where X^ is the x drop position and is the z drop position. 

The integration of a droplet trajectory begins at a point suf- 
ficiently forward of the airfoil where the flow is unaffected by 
the airfoil and continues until either the droplet impacts the 
airfoil or passes over or under the airfoil. 


4.1 OVERALL COLLECTION EFFICIENCY 

The overall collection efficiency for a given drop diameter 
is the ratio of drops that strike the airfoil to those that would 
have hit the airfoil had the airflow not affected the drop 
trajectories. This ratio can be expressed as E ■ (Z -Z )/H where 

n ii 

H is the projected height of the airfoil onto the initial drop 
positions and Z„ and Z are the initial drop positions for 
respectively the highest and lowest drops to impact the airfoil 
(see Figure 4). 

Some examples of trajectories in the vicinity of the wing 
are shown in Figures 5/ 6, and 7 together with the overall 
collection efficiency of each drop size. The respective drop 
diameters for these figures are .1 mm# 1 mmr and 3.1 mm. 

The overall collection efficiency by rainfall rate is 
obtained by summing the product of drop size collection effi- 
ciency and percent volume contribution for all size drops. The 
summation is done for drop diameters from .5mm to 8mm in .5mm 
increments. 

16 

CE • I E.M(D). 
i-1 ^ ^ 

where E^ and N(D)^ are respectively the overall drop collection 
efficiency and percent rain volume of the ith size drop. Table 1 
summarizes the results for a symmetric airfoil at 0* angle of 
attack with a 65 m/sec velocity. The overall collection effi- 
ciencies are greater than 95% for all rainfall rates of Interest. 
Similar results would be anticipated for other airfoils and for 
the fuselage. Large angles of attack may decrease somewhat the 
overall collection efficiencies. 
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Figure 4, Tangential Drop Trajectories for Potential Airflow About an Airfoil 
Drop trajectories begin well ahead of the airfoil to ensure that 
initially their paths are unaffected by the airflow. 
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Piguire 6. Trajectories of Iran Diameter Raindrops in Potential Flow About an Airfoil. 












TABLE 1 

OVERALL IMPINGEMENT EFFICIENCY BY RAINFALL RATE 


Drop 

Diameter (mm) 


Collection 
Efficiency (%) (E) 


0.5 

78.4 

1.0 

87.1 

1.5 

95.2 

2.0 

96.4 

2.5 

97.7 

3.0 

99.0 

3.5 

99.5 

4.0 

99.5 

4.5 

99.5 

5.5 

100.0 

6.0 

100.0 

6.5 

100.0 

7.0 

100.0 

7.5 

100.0 

8.0 

100.0 



CE » 



(%Vol)^ 


% Vol is volume contribution 
of ith size drops 


Rainfall Rate (mm/hr) 2000 1000 

E (%) 99+ 98.7 


500 

98.1 


200 

97.4 


100 

96.7 


25 

95.7 
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4.2 LOCAL COLLECTION EFFICIENCY 


ORIGINAI. PACE !3 


Determining where a raindrop has impacted the wing is 
crucial to calculating the local collection efficiency at various 
segments of the airfoil as well as calculating other impact 
parameters. Local collection efficiency is required in assessing 
the water film that may develop on the wing due to heavy rain in 
addition to assessing the roughness that may develop in the film. 
The local collection efficiency^ i, is the initial vertical 
distance between successive drop trajectories « dt, divided by 
the distance along the wing between successive impacts f ds (see 
Figure 8) . Using the droplet trajectory program* local collec- 
tion efficiencies were calculated for each droplet diameter at 
various stations around the airfoil. Local collection efficien- 
cies by rainfall rate were then established at each st^/tion by 
weighting the percentage volume contribution of each droplet 
diameter by the local collection efficiency for that diameter and 
summing over drop diameter. 



Figure 8. Calculation of Local Collection Efficiency. 
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SBCTION 5 

MONBMTUH OP RAINDROPS 

Raindrops striking an aircraft loss noasntua to ths 
aircraft, thus changing ths vslocity of ths aircraft, 'ttis ver- 
tical conponsnt of ths raindrop vslocity iaparts downward wwsn— 
tun to ths aircraft# which tends to sake it sink. Ths raindrops 
striking head-on slow ths aircraft because nonsntum is lost in 
accelerating ths water droplets to ths aircraft vslocity. Wis 
aaount of noiasntua imparted to an aircraft by striking a raindrop 
depends on ths reflection angle of the raindrops. A raindrop 
striking the aircraft surface acutely iaparts less aoMntua than 
if it strikes at a larger reflection angle. To estiaate the 
iapacted aomentun for a large aircraft such as a Boeing 747, the 
following assuaptions are aades 

(a) All rain iapinging on the aircraft accelerates to the 
velocity of the aircraft (inelastic collision) 

(b) The aircraft is in straight and horisontal flight at 
sero degrees angle of attack 

(c) The flaps are not deployed 

(d) The aircraft goes froa the no-rain to rain situation 
instantaneously . 

The assuaption that the rain strikes the aircraft directly 
and takes on its velocity is an underestiaate for drops striking 
head on. Any reflection for these drops will be at an angle 
greater than 90oi thus, the iapact force froa the partially 
elastic collision of these drops is actually greater than 
accounted for in the inelastic calculations. Boas drops striking 
the wing or fuselage will suffer reflection at an acute angle. 

The force iaparted by these latter drops is less than that 
accounted for in the aoaentua calculations. The assertion that 
the iapacts are inelastic is supported by the Lucey (1972) work 
on the Crush of a Point Detonating Blexnt. 7?he eleaent was 
mounted on a rocket sled and passed through a heavy rainfield. 
The naasured crush of the element was related to work energy and 
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wa« ooapared to the work that would have been dona by both 
*^**^^® ••'<5 inalaatic oolliaiont. Tha eoapariaon ahowad tha 
maasurad crush could ba bast axplainad by inalaatic collisions. 

Tha calculation of rain-inducad monantum panaltias doas not 
consider a landing wing configuration using flaps; this wnuld 
P*^****'^ • larger inpact cross-section than doas a wing in hori- 
zontal flight. Likewise/ a non-horizontal aircraft attitude 
would generally present a larger inpact cross-section. A greater 
nomantun penalty would result in both cases. 

The above assunptions used in calculating nonantun 
panaltias nay ovarastinata tha affect in sona areas and 
underest inata it in other areas. Nonetheless/ we believe they 
exhibit sufficient realise for an initial aerodynanic assassnant. 

To estinate tha horizontal nonantun inpartad to tha 
aircraft/ visualize the process as that of a water jet whose 
discharge is equivalent to the rainfall interception rata. Tha 
water jet hits tha aircraft at tha fraastraan spaed and loses all 
its nonantun to tha aircraft. What is tha force tha aircraft 
nust apply to ovarcone tha jet? 

In tha horizontal x-diraction 




where 


H is tha nass interception rata of water par second/ «nd 
^xPREE horizontal conponant of tha aircraft velocity rala- 

tiva to tha air nass. With tha assunption that tha final velo- 
city of tha water after inpact is zaro/ Bq. (7) datamines t!ia 
horizontal force exerted on tha aircraft by tha rain. 

Tha vertical force of rain inpact on tha aircraft is con- 

lass since tha taminal velocity of heavy rain in still 
air fron aq. (2) is about 9 n/sac. 
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In the vertical z-direction 


F. 


" <"o ♦ '',pree' 


( 8 ) 


where 

is the vertical velocity o£ the drops , and 
V^Pree is the vertical component o£ the aircraft velocity 
relative to the air mass. 


The mass interception rate, A, is given by 

M • iV^RlApPjjCE 

where 

v5r » Impact velocity vector of water drops 

A * Projected area o£ the aircraft in the direction of 
P 

drop impacts, Wr 

p « liquid water content of the air as a function of 

<w 

rainfall rate 

CE • overall collection efficiency of the aircraft as a 
function of rainfall rate and angle of attack. 

The impact velocity vector of the water drops is derived 

from summing the Aircraft Velocity Vector (V ), the Wind Vector 

AC 

(V^jr) , and Raindrop Terminal Velocity Vector (W^) as 


Wr - V - V^iR * W^- 


The projected area A^ for an aircraft varies considerably 
as the view angle changes from horizontal. An accurate deter- 
mination of the presented view of a specific aircraft by view 
angles is beyond the scope of this report. A flat surface 
approximation can be used to obtain the projected area as a func- 
tion of view angle, given the frontal and top view surface areas. 


Th« approximation, is 

Ap • A^sin(if) + ApSin (90o-t) 

where 

• Top view surface area 

A„ ■ Frontal view surface area 
F 

t » Angle between fuselage reference line and 

Evaluation of the total force exerted on a 747 aircraft and 
the X and z components of this force for horizontal flight at 65 
m/sec (approximately 125 knots) and 0* angle of attack is given in 
Table 2 by rainfall rate. The additional force needed to balance 
the momentum penalty can be compared with the thrust produced by 
the aircraft engines. For a 747 aircraft the maximum engine 
thrust is on the order of 800,000 newtons (180,000 lbs). Thus at 
a rainfall rate of 100 nm/hour only 0.4 percent of maximum thrust 
is needed to counteract the rain momentum while at a 2000 rom/hour 
rate, 9 percent of maximum thrust is required. 

If no additional thrust were applied and the other forces 
on the aircraft remained constant then the rain momentum would 
extract speed from the aircraft. The resulting deceleration 
equals the momentum force divided by the mass of the aircraft. 

The deceleration for a 180,000 kg aircraft varies from .04 
knots/sec for a 100 mm/hour rain to 0.2 )cnots/sec for a 500 
mra/hour rain to 0.75 knots/sec for a 2000 mm/hour rain. If the 
aircraft were in the heavy rain environment for a 20 second 
period, the approximate resulting loss of airspeed would be 0.8 
knots for the 100 mm/hour rate, 4 knots for the 500 mn/hour rate 
and 15 knots for the 2000 mra/hour rate. These calculations of 
airspeed loss could be an underestimation for a landing con- 
figuration with high lift devices extended (which increase the 
water catch rate) or for a landing at a higher approach velocity 
(larger momentum loss) or when applied to an aircraft executing a 
go-around maneuver (increased water catch rate). Nevertheless it 
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appears that significant momentum penalties result for transport j 

class aircraft for rainfall rates of 500 mm/hour or greater. For 
rates less than 500 rom/hour » typical of roost accident scenarios , 
the momentum penalty may be a contributing factor but of itself 
would not be expected to present severe problems. These momentum 
penalties are further analyzed in a subsequent section using a 
landing simulation program to assess their effect on aircraft 
performance. 


TABLE 2 


FORCE 

EXERTED 

ON AIRCRAFT DUE TO 

MOMENTUM OF 

DROPS 

RAINRATE 

‘’iw 

W 

O 

F 

P , 

A 

F 

z 

(mm/hr) 

3 

gm/m'* 

m/sec 

Newtons 

Newtons 

Newtons 

100 

3.23 

8.42 

3.60x10^ 

3.57x10^ 

4.57x10^ 

200 

6.23 

8.96 

7.13x10^ 

7.06x10^ 

9.67x10^ 

500 

15.31 

9.14 

1.82x10^ 

1.80x10^ 

2.62x10^ 

1000 

30.18 

9.30 

3.58x10^ 

3.54x10^ 

5.17x10^ 

2000 

59.74 

9.45 

7.09x10^ 

7.01x10^ 

1.02x10^ 



B-747 

Aircraft 





\c - 

65 m/sec 





\ " 

1131 m^ 





""f " 

119 m^ 




CE - 1.0 
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SECTION 6 
WATER FILM 


In very heavy rain, a water film forms on the upper sur*> 
faces of the wings, fuselage, and tail. The film may profoundly 
affect aircraft performance not only by increasing aircraft mass, 
but also by increasing total aircraft drag and decreasing lift. 
For the purposes of aerodynamic penalties, it is necessary to 
calculate the film's mean thickness. In order to arrive at the 
film thickness, a number of assumptions were made. 

Because it is relatively thin and thus moves considerably 
slower than the airflow above it, the film is assumed to be lami- 
nar. Its primary motivation is due to the airflow. At the 
interface viscous stress is matched between air and water. The 
resulting film is a balance between water runback and rain water 
reception. We assume the rain water reception. is not affected by 
droplets shedding from the drop impact crowns. The shed droplets 
are probably returned to the film by boundary layer entrainment. 

The local rate of mass reception is required for calcu- 
lating the film thickness. The reception rate is dependent on 
the liquid water content of the air, the free stream 

airspeed, th® local collection efficiency 0. For a 

given location the reception rate per unit area is given by 


M « 


free* 


A numerical model for the water film was developed to 
caJculate film thickness by airfoil location and rainfall rate. 
The model geometry calculates the film flow at a number of sta- 
tions along the wing as shown in Figure 9. At each station, the 
film thickness results from water remaining after the gains and 
losses due incoming rain and film flow are considered. After a 
time, an equilibrium thickness results at each station. A more 
complicated model that considers pressure gradient, gravity, and 



ORIGINAL FA3E IS 

OF POOR QUALITY 






other effects U not werr.ntea In a first look at heavy rain 
effects on aircraft. Refinements to the model may be desirable 

in the future however. 

For simplicity, the model considers the wing to be a flat 
plate and the direction of the arriving rain varies from nearly 
normal at the leading edge to nearly tangential depending on 
angle of attack at the trailing edge. Slower airflow near the 
stagnation point is not Incorporated, for this reason the calcu 
lated film thicknesses near the leading edge are probably too 


thin. 

The model is two-dimensional and consists of computations 
at 100 stations shown in Figure 10. A linear velocity profile is 
assumed between specified velocities at the film's top and bot- 
tom. At the bottom, the film velocity is zero. At the top, the 
velocity is based on the viscous stress due to the 
above the water film. Following Hartley and Murgatroyd (19«4), 
the shear stress at the surface of the film is 


, . f.pV.Vl 


<») 


where f is the air friction coefficient, p is the air density 
and Vj is average air velocity over the wing. The velocity in 

the film at height * is 

(10) 

u ■ xz/m 


where u is the fluid viscosity. 

The flow was calculated by use of the u momentum equation. 


du 

dt 



( 11 ) 


and the continuity equation. 


aw . 
az 


au 

az 


( 12 ) 
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Figure 10. Hater Film Thickness Model (2). 
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The calculation begins at the second station. Vertical 
velocities r Wj are obtained from eg. (12) and u velocity tenden~ 
cies from (11) beginning at one point below the film top and pro- 
ceeding downward to one station above the bottom. This procedure 
is repeated sequentially through the 99th station and its comple- 




tion constitutes one time step. Using a sufficiently small time 
step to ensure computational stability# the u velocity at time t 


•f At is calculated from 


2 

U(t+At) - U(t-At) ■ 2At -U*'* -W^“ 2 

3x 9z az^ 

For the first time step a forward difference is used. Two 
neighboring stations malce up a volume into which and out of which 
there is mass flux due to the arriving rain and film flow. 

The mass balance at a station# i# may be expressed as 


dm^ 

Tt 


p z«h 
^ / udz 

^ Z“0 
X-iAX 
( 1 ) 


p z«h 
^ / udz + 
^ z»o 

X«u+1) AX 
(2) 


ami 

at 

(3) 


(14) 


where terras 1 and 2 represent respectively mass flux into 
the ith station (x»iAx) of Area A and out of the station 
(x>(i-t'l) AX) due to film flow. Term 3 represents mass flux# into 
the ith station due to rain. 


A change of mass is calculated from Eg. (14). 
A new film thiclcness# h^# results from 


m^ 


dm^ 

• 

P A 


At 


(15) 


The procedure is repeated for each station from first to last. 
The calculations were performed from an inicial time t«0 with 
an initial m^^ corresponding to an initial thickness of 0.5 nm 
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until steady state was established. Other initial values of 
did not significantly change the steady state value. 

The resulting average film thicknesses for the top of the 
airfoil derived by the model by rainfall rate are summarized in 
Table 3. The calculations were made for a symmetric airfoil of 
chord length 10 m at 0* angle of attack and 65 m/sec (126 knots) 
free stream airspeed. The value shown in Table 3 it* an average 
representation of the film thickness throughout the airfoil 
segments. 

The film thickness for the fuselage was estimated from 
calculations for an airfoil. Because of increased surface area 
of the fuselage in relation to its projected area, the droplet 
impact density on the fuselage is approximately 2/3 that on the 
airfoil. Thus, the mass impact rate per unit area on fuselage is 
only 2/3 that of the airfoil. The fuselage film thickness at a 
given rainfall rate is estimated as the airfoil film thickness 
that corresponds to 2/3 the given rainfall rate. 

TABLE 3 

AVERAGE FILM THICKNESS FOR A SYMMETRIC AIRFOIL AND 
FUSELAGE AT 0* ANGLE OF ATTACK, 10m CHORD 


Rainfall Rate 
(mm/m) 

Calculated 
Thickness 
Airfoil (mm) 

Estimated 
Thickness 
Fuselage (mm) 

100 

<N 

• 

O 

VI 

<0.2 

200 

<0.5 

<0.2 

500 

0.8 

O.fi 

1000 

1.0 

0.9 

2000 

1.3 

1.1 
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SECTION 7 
WEIGHT PENALTY 

The weight penalty is due to a water film on the surface of 
the aircraft. Using a 747 as an example with the average film 
thicknesses calculated in Table 3# the weight penalties by rain- 
fall rate are shown in Table 4. The values are for a film 
existing only on the top half of the wing and fuselage surfaces. 
Film thickness was not calculated for the lower surface, below 
the stagnation point. If a film exists there, it would be 
expected to be thinner because of a decreased water impact rate. 
Even assuming the film on the underside was equally thick it 
would only double the weight shown in Table 4. For a 747 
aircraft whose landing weight is on the order of 180,000 kg, the 
maximum weight penalties cannot be much larger than 1% of landing 
weight. A 1% added weight has a negligible effect on aircraft 
landings. 


TABLE 4 

WATER MASS ON ALL AIRCRAFT SURFACES FOR VARIOUS RAINFALL RATES 


Weight Penalties 


Aircraft ■ B747, 

Velocity » 126 

kts, a - 0* 

Rainfall Rate (ma/hx) 

Film weight, 
top (kg) 

Film weight, 
top and bottom 

100 

— 

— 

200 

7.33 « 10^ 

1.47 * 10^ 

500 

1.17 >c 10^ 

2.34 * 10^ 

1000 

1.47 K 10^ 

2.94 X 10^ 

2000 

1.94 * 10^ 

3.88 X 10^ 
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SECTION 8 
AIRFOIL ROUGHNESS 

An airfoil or fuselage in heavy rain may be roughened in at 
least three ways: (a) waves develop in a water film clinging to 

the wings and fuselage, (b) in the absence of a liquid film, glo- 
bules dot the wing surface and are blown back by wind stress, (e) 
drops impacting a water film disturb its surface. 

This study has analyzed two of the roughness sources. They 
are the roughening of a water film on an airfoil's top surface due 
to d«:'op impacts and the roughness due to waves in the water film. 
Future modeling can consider the possibility of blow back of 
liquid water globules. 

8.1 ROUGHNESS DUE TO IMPACT 

We have investigated the impact of a raindrop on a thin 
film in order to evaluate its contribution to roughening an 
airfoil. The assessment of raindrop splashes is based on the 
work of Machlin and Netaxas (1976) %rho investigated the craters 
formed when a water drop impacts a thin water film. 

Unfortunately their investigation was conducted at drop veloci- 
ties slower than the impact velocity of raindrops on an aircraft 
wing. Thus extrapolation of their results to higher velocities 
was required. 

The raindrop splash model assumes that a drop hitting a 
thin film forms a cylindrical crater. All the water originally 
residing within the crater is assumed to go into the crater 
crown. Virtually no wave swell results from the impact. An 
energy balance equation relating the non-dimensional crown height 
to crown radius in terms of the inertial, gravitational (Froude 
Number) and surface tension (Weber Number) is given by: 


(Macklin and Mataxaa, 1976) 


* 

(H ) 


1 ♦ - 3(Wjj)"^D*R* 


( 16 ) 


whara 


w. 




* 


D 

* 

H 


« 


R 


C 


* 

t 



0 




h/Rj 

H/Rj 

V»d 


We bar Numbar (ratio of inart ial to 
aurfaca tanaion forea) 

Prouda Nuabar (ratio of inartial to 
gravitational forca) 

Dinenaionlaaa dapth of liquid film 

Oinanaionlaaa crown haight 

Oinanaionlaaa crown radiua 


Vfraa^/^d Oinanaionlaaa tina 

drop radiua 
aurfaca tanaion 


Thia aquation doaa not taka into account viacoua forcaa 
(which contribute laaa than 7 parcant) and other influancaa not 


accounted for in thia ideal ixad nodal. In oonparing thia nodal 
with axparinantal data Macklin and Nataxaa found the nodal to 
account for 40 to 60 parcant of the energy diaaipation in ahallow 
filn inpacta. The unaccounted energy ia thought to pradoninataly 
conaiat of aurfaca energy. The aolution of Bq. (16) raquiraa the 
uaa of an ^irical ralationahip between dinenaionlaaa crown 
haighta* H , and crown radiua^ R *. Pigura 11 ahowa the ratio of 

W 

* * 

B different Nabar nunbara. It waa nacaaaary to axtrapo* 

lata the curve fron Pigura 11 to nuch higher Mabar nunbara for 

ai^lication to dropa inpacting on an aircraft. An extrapolated 

ratio of 2.9 waa uaad for a Nabar Nunbar of approxlnataly 

5 * 10^. Uaa of auch a aavaraly extrapolated value waa unde- 

airabla but nacaaaary. Sona conaolation can be gained by a aan- 

aitivity analyaia uaing other ratioa of H/R^. Uaing ratioa of 

c 
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Linear Regression Fit 
Macklln and Netaxas Data 



Ratio of Cro%m Height, H, to Cavity Rdaius at Maximum Cro%m Height 
as a Function of Weber Number using Macklin and Metaxas (1976) data 








H/R ■ 1 and H/R^-5 doaa not change the aasantial results o£ the 
analysis. Using the extrapolated ratio, H is theoretieany 
obtained fron the solution of Bq. (16). The theoretical crown 
height is then adjusted by Figure 12 to account for the energy 
dissipated that is unexplained by the nodel. 

The dinensionless time to reach maximum crown height is 
obtained from Figure 13 as a function of experimental and theo- 
retical values of H . Figures 11, 12, and 13 each required 
g^v 0 pe extrapolation of the experimental data for the situation 
of rain impacting on an airplane. 

The values derived for maximum crown height, radius, and 
time to reach maximum height are the basic parameters for 
deriving the aerodynamic roughness height due to drop impacts on 
an airfoil. Table 5 shows these parameters for a drop impact 
velocity representative of a transport aircraft landing. 

TABLE S 

DROP IMPACT PARAMETERS BY DROP SIZE 
V«65M/sec 


Drop 

Radius 

(am) 


'n 

H/*' 

* c emp 

# 

®theor 

**emp 

* 

*e 

0.5 

2.98>10^ 

8.61*10® 

2.79 

113.94 

13.2 

41.5 

1.0 

5.79*10^ 

4.31*10® 

2.98 

169.48 

14.6 

43.4 

2.0 

1.16*10® 

2.15*10® 

3.25 

250.35 

15.9 

45.0 

4.0 

2.31*10® 

1.08*10® 

3.53 

368.71 

16.8 

46.0 


Evaluating the roughness of a film due to drop i^acts 
requires Icnowing two parameters. One parameter is the mean 
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Experimental Results III 



Figure 13. 


The Dimensionless Time to Reach Maximum Crown Height 

t* as a Function of the Theoretical and 

crown 

Experimental Values of H* for Shallow Liquid 
Splashing. (After Macklin and Metaxas, 197 d ) . 
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height of a drop impact crown throughout its life cycle from 
formation to dissipation. The other parameter is the mean 
separation distance between drop impacts. Both parameters can be 
calculated from the results in Table 5 and local collection 
efficiencies. 

The mean separation distance is derived from the rate of 
drop impacts and the time over which a drop impact persists. 

The rate of drop impacts of diameter D at a given location on the 
airfoil is the product of drops per unit volume times the 
freestream air speed times the local collection efficiency. The 
time periods At» over which a drop impact crater exists is con- 
servatively approximated as twice the time required for a crater 
to reach maximum crown height, that is: 


It • 


(17) 


The number of drop impact crowns of diameter D occurring at any 
given time is 

N(D)Vfree^t® » N(D)t*$D (18) 

where N(D) is number of drops of diameter D per unit volume in 
free air. Using the Marshall -Palmer distribution to represent 
N(D) and integrating the above over all D gives the total number 
of drop impact craters at a given time, 

* 

D*** ,n * 

N - / N e’^^t^fiDdD • . (19) 

Dio o c *2 


In integrating eq. (19) both t* and 6 are assumed constant as 
they are slowly varying functions of 0. 

The mean separation distance is given as the inverse 
square root of the number of drops involved in an impact at 
any instance of time per unit area. That is: 




For a drop of diameter D, the average drop impact crown 
height over its lifetime is conservatively approximated as one 
half of its maximum crown height. That is: 


k(D) . (21) 

The average height, R, for a population of drop impact cra- 
ters of varying diameters is the population weighted mean 


/ k(D)N(D)dD 

D«0 

D* ^ 

/ N(D)dD 
D-0 


( 22 ) 


The evaluation of eq. (22) again assumes t^ and $ are constant 
and that N(D) is represented by a /larshall -Palmer distribution. 
Thus, 


H 0 
2 * 


(23) 


where k is considered to be the height of the roughening elements 
produced by drop Impacts. 

8.1.1 Sandqrain Roughness 

The roughness height is estimated from Eqs. (20) 
and (23) by finding the equivalent sand grain roughness (Dlrling, 
1973). We first found the ratio A^/A^, where is the projected 
area of drop crowns in the freestream velocity direction and A 
is the windward surface area of the element as seen by the flow. 
We estimate this ratio to be about 0.64 for cylindrical drop 
impact crowns. Prom Dirling, A is correlated in terms of the 
spacing parameter, D/k, as 

A - (D/k)(Ap/A^)“^/^ (24) 


Considering, for simplicity, an average local 
collection efficiency of 6 ■ .17 for all rainfall rates, and 
using the value of t^ corresponding to the mean drop size for 
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a given rainfall rate» D was calculated from (20) for a free 
stream velocity of 65m/sec and is given in Table 6. 

TABLE 6 

AVERAGE SPACING BETWEEN RAINDROPS BY RAINFALL RATE 





Rainfall 

rate (mm/hr) 

D (cm) 


100 

20.16 


200 

17.43 


500 

14.38 


1000 

12.43 


2000 

10.74 

Using 

the same assumption. 

the average geometric 

height of drop impacts was calculated for 

different rainfall 

rates using Eq. (23) 

as shown in Table 7* 



TABLE 7 


AVERAGE GEOMETRIC 

HEIGHT OF DROP IMPACTS BY RAINFALL RATE 

Rainfall 

rate (mm/hr) 

Height k(mm) 


100 

4.8 


200 

5.6 


500 

6.7 


1000 

7.8 


2000 

9.0 


From Dirling the correlation equations for 
sandgrain roughness developed from experimental data are: 
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* 0.0164A 


3,7 


/ A< 4.93 


( 25 ) 


13 9 A 


-1.9 


r A> 4.93 


where ■ equivalent sandgrain roughness. 

The average sand grain roughnesses by rainfall rate 
derived by eqs. (24) and (25) are given in Table 8. 

TABLE 8 

SAND GRAIN ROUGHNESSES BY RAINFALL RATE 



8.1.2 Drag Increase Due to Impact Roughness 


To make estimates of the drag increase due to these 
values of sandgrain roughness, results from Young (1965) for fixed 
roughness elements with turbulent flow over a flat plate were 
utilized. The mean friction coefficient, , for smooth wall 
flow is given by: ^ 


Cp ■ 0.088/(log Re-1.5) ^ 

S 

In roughened flow the mean friction coefficient is: 
Cp - (1.89 + 1.62 log L/k.)"2.5 


(26) 


(27) 


where L is the mean aerodynamic chord or fuselage length. 

For a 747 wing, an appropriate Reynolds number in the landing 
configuration is Re - 3.23 x lo^ while for the fuselage it is 
2.74 X 10®. 
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Using values of 8.3 meters for the Mean Aerodynamic 
Chord of a B747 and 70 meters as the fusexage length# the solu- 
tion of Eqs. (26) and (27) are shown in Table 9. 


TABLE 9 

MEAN FRICTION COEFFICIENT FOR SMOOTH AND ROUGHENED 
AIRFOIL AND FUSELAGE 


L * 

8.3m Airfoil 

L » 

70m Fuselage 


Rain Rate 

Cp 

's 

Co 

^R 

^F 

^S 

Cp 

*R 

mm/hr 

Airfoil 

Airfoil 

Fuselage 

Fuselage 

100 

.0024 

.0036 

.0018 

.0025 

200 

.0024 

.0042 

.0018 

.0028 

500 

.0024 

.0051 

.0018 

.0033 


.0024 

.0059 

.0018 

.0038 

2000 

.0024 

.0069 

.0018 

.0043 


An estimate of the influence of an increase in 
friction drag on total drag for the 747 aircraft in the landing 
configuration can be approximated as follows. For a 747 aircraft 
with 30 degrees flaps descending a glide slope at 2 degrees angle 
of attack:# the basic drag coefficient including landing gear is 
approximately C^ « .15. The basic drag coefficient can be 
decomposed into Sontributions from airfoil and fuselage friction 
coefficients as: 


where 



^S ^S 


+ c 


D etc 


(28) 


fus 


'D etc 


total surface area of the fuselage 

wing area (upper surface) 

all other factors that contribute to , 
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Thtt factor, 2, in tha first term on th« right hand aide accounts 
for friction drag on both upper and lower wing surfaces, k 
representative value for the ratio of fuselage to upper wing 
surface area for a 747 aircraft is 


A 

S 


fus 


3.4 . 


For an aircraft whose upper surfaces are roughened 

by rain, the change in C due to increased friction drag is< 

o 

.fus 


AC« « 


pair , pair A^^^pfus pfus 


(29) 


Using the appropriate values from Table 9 in Bq. (26) the percent 

changes in C for a 747 aircraft in the landing configuration 
o 

are shown in Table 10. 


TABLE 10 

INCREASE IN TOTAL DRAG DUB TO INCREASED NING AND FUSELAGE 
FRICTION DRAG. (747 AIRCRAFT LANDING CONFIGURATION) 


Rainfall Rate 
(nn/hr) 

D D- 
o 

100 

1.6% 

200 

2.3% 

500 

3.5% 

1000 

4.6% 

2000 

5.9% 


The increases in drag coefficient shown in Table 10 
are small but significant. If both the upper and lower surfaces 
of fuselage and airfoil were roughened to the same extent, the 
values shown in Table 10 would increase by a factor of 2. 
Similarly, if the above analysis were performed on a smaller 
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aircraft r such as a 727 or DCS a larger percentage increase in 
friction drag would be derived from eqs. (26) and (27). 

To establish the validity of applying eq. (26) and 
( 27 ), which refer to flat plate turbulent flow conditions to an 
airfoil with extended flaps, a comparison of theory to experi- 
mental wind tunnel measurements by Ljungstroem (1972) was made. 
Ljungstroem measured both lift and drag increments on a 
2-dimensional wing section with and without high lift devices. 

The wing chord length was « 65 cm. Figure 14 shows and 
curves for a clean wing and a roughened wing (95% and 100% 
coverage) with various sand grain roughness elements. Using 
eqs. (26) and (27), the appropriate Reynolds number, L <■ 65cm 
and )c « 0.1mm, k • 0.5mm, the smooth and roughened friction 
coefficients were calculated as C_ » .0036, C_ « .0052 (k • 

.1mm) and C « .0076 (k. ■ .5mm). The values are shown in 
^R * 

Figure 14 after converting to drag coefficient. Note that 
Bq. (26) gives a good approximation to the smooth airfoil drag 
coefficient but (27) grossly underestimates the roughness 
penalty by a factor of 3 or more. Since the wind tunnel measure- 
ments include pressure drag as well as friction drag it is 
apparent that additional drag penalties result from roughness 
above those calculated by Eq. (27). Similar results «#ere 
obtained from a comparison of Ljungstroem measurements on a high 
lift airfoil with 25-degree slats and 20 degrees flap. Figure 15 
shows the measured vs at low angles of attack as well as 
the friction drag calculations using Eqs. (26) and (27). The 
experimentally derived increase in between the clean and fully 
roughened airfoil is approximately • 0.025 for k^ « 0.5mm and 
» 0.009 for O.lima roughness. Calculated differences in skin 

friction coefficient are AC_ « .0040 for k - 0.5mm and AC > 

r a F 

.0017 for k • 0.1mm. For this wing configuration the skin fric- 

D 

tion calculation underestimates the roughness inJuced drag 
increase by a factor of 5 to 6. Thus in applying the increase in 
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(Abbot & Ooenhoiff Ref 3) 


Figure 14. Cj^ and for Configuration W. 

Wing Section NACA 65 A-215. (After Ljungstroem; 
1972) . 
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Figure 15. C^j vs for Configuration WFS 

Wing Section with Slat and Training Edge 
Flap. (After Ljungstroem; 1972) . 
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friction drag calculated by Eqa* (26) and (27) aa the aole 
contributor to increased drag coefficient and underestimation by 
a factor of 3 to 6 appears li)cely. Consequently the percentage 
drag penalties due to roughness shown in Table 10 should bs 
increased by a compensating factor - perhaps on the order of 3* 

8,1,3 Summary of Impact Drag Penalties 

The above results on drag penalties that arise from 
the roughness of an airfoil and fuselage due to rain drop impacts 
indicates aerodynamic penalties may exist. Numerous assumptions 
and experimental data extrapolation however were required to 
derive the overall drag penalty for a 747 aircraft. Most criti- 
cal of these assumptions relate to the extrapolation of droplet 
crown height and the time to reach maximum crown height for Weber 
numbers far beyond the range of experimental measurements. Also 
the assumption of fixed roughness elements, no Interaction of the 
airflow with the crown formation, and no shattering of the water 
drops prior to impact, may be important. Nonetheless the derived 
roughness heights when evaluated in comparison with experisMn- 
tally derived drag penalties appear especially significant. 
Furthermore since the roughness due to the waviness of the water 
film has not yet been taken into account, a further contribution 
to the drag penalty may result. Finally, if the same analysis 
were performed on a smaller size transport aircraft there is 
reason to believe the drag penalty may Increase. 

8.2 ROUGHNESS DUE TO FILM WAVINESS 

Surface waves are observed on water film surfaces under 
wind stress (Kapitsa. 1964). They have been shown to have a com- 
parable effect on airflow as sand roughened walls. Wurs (1978) 
conducted an experimental program to study the interference of a 
wavy liquid film with a turbulent gas boundary layer. Using a 
two-phased wind tunnel of .72 meters in length, and appropriate 
instrumentation. Wurz measured various properties of the water 
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film and its interface with the air boundary layer for a range of 
Mach numbers from subsonic to supersonic* Of interest to our 
application for rain roughness on an airfoil are the results at 
low subsonic Mach numbers. Tests conducted in the Mach number 
range of Mach ■ *18 to Mach « *6 produced film thickness in the 
range of .014 mm to .2 mm. The water flow rates associated with 
these film thicknesses are an order of magnitude less than that 
produced by an aircraft penetrating even a 300 mm/hr rainfall 
rate at landing speed. Thus# the results of the Wurz experimen- 
tal work# even though it covers sufficiently low Mach numbers# 
must be extrapolated to thicker films than occurred in the 
experimental tests. Wurz's results however# indicate an increase 
in friction coefficient as the film thickness increases. A rela- 
tionship was also established between the mean film thickness and 
equivalent sand grain roughness as derived from the roughness 
Reynolds number; this ratio varies from between 1.25 and 2.0 in 
the subsonic test range. Using a ratio of 1.5# an equivalent 
sand grain roughness was derived for the airfoil and fuselage 
film thicknesses associated with 100 mm to 2000 mm rainfall 
rates. Table 11 shows these results. 



TABLE 11 

EQUIVALENT SAND GRAIN ROUGHNESS OP NAVY WATER FILM 


Rain Rate Equivalent Sand Crain Roughness k^ 

(mm/hr) Airfoil (mm) Fuselage (mm) 


I 




Again using aq. (27) r astiaatad tha incraaaad skin 
friction coafficiant. Tabla 12 shows thasa valuas for tha 
airfoil and fusalaga varsus rainfall rata. 

7ABLB 12 

INCREASE IN SKIN FRICTION DUE TO FIUI NAVINE8S 


Rain Rata 
(m/hr) 


'R 

Airfoil 


'R 

Fusalaga 


200 

.0048 

.0032 

500 

.0053 

.0032 

1000 

.0055 

.0036 

2000 

.0059 

c0038 


Using aq. (29) thasa incraasas in friction drag froai Tabla 
12 wara convartad to par cant incraasa on drag coafficiant for 
wavas on tha uppar aurfaca of tha wing and fusalaga. Tabla 13 
shows thasa rasults. 
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TABLE 13 

INCREASE IN DRAG COEFFICIENT DUE TO FILM WAVINESS 


Rainfall Rate 
( mm/hr ) 

(») 

100 

2.1 

200 

3.2 

500 

3.8 

1000 

4.2 

2000 

4.6 


Drag increases in the range of two to five percent are 
derived for the various rainfall rates. This drag increase is 
based solely upon an increased friction coefficient derived from 
Eqs. (26), (27), and (29). If» however, the drag coefficient 
penalty was derived using the sand grain roughness calculation of 
Table 11, and the experimental measurements of Ljungstroem, then 
the drag penalties would increase by a factor of from two to six 
over chose of Table 13. Thus there is reason to suspect that 
Table 13 may seriously underestimate the actual drag penalties 
associated with the waviness of the water film. 

in comparing the derived drag penalty due to film waviness 
(Table 13) with that derived for drop impact cratering, (Table 
10) it is seen that the impact cratering penalty is approximately 
the same as that associated with film waviness. However, the 
experimental data upon which the waviness calculations were based 
require less severe data extrapolation and thereby more con- 
fidence is placed in these values. In combining the drag estima- 
tes due to waviness with those due to drop impact cratering and 
ta)cing into account the possible underestimating of these 
penalties when compared to experimental measurements, a best 
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overall estimate of total drag penalty for a transport aircraft 
in the landing configuration was made. We estimate this penalty 
for 100 to 2000 mm/hr rainfall rates to be in the range of ten 
percent at the lower rainfall rates up to 20 percent or more for 
the incredible rates. 

8.3 REDUCTION OP LIFT DUE TO ROUGHNESS 

Distributed roughness on the upper surface of an airfoil 
also affects lift. Brumby (1979) summarized the results of 23 
experimental investigations concerning the effect of toughness 
on lift coefficient and stall angle. Included in Brumby's analy- 
sis was the data of Ljungstroem for airfoils with and without high 
lift devices (see Figure 4). For a fully rough airfoil with 
"large* roughness elements, the lift coefficient was found to 
decrease at all angles of attack. For "medium" roughness a 
decrease in the lift coefficient primarily occurs at high angles 
of attack. The decrease in maximum lift coefficient (CL ) is 
related to the ratio of roughness element height to wing"chord. 
Decreases in CL^^^ as high as 30-40 percent occur for "large" 
roughness elements. A decrease in stall angle also accompanies a 
decrease in A stall angle decrease of 4 to 5 degrees is 

appropriate to a decrease in CL^^^^ of 30 to 40% while one to 
three-degree decreases in stall angle is appropriate for lift 
losses of 5 to 15%. An increase in stall speed also occurs. 

Stall speed increases of 10 to 20 knots accompany a loss in maxi- 
mum lift of 10 to 30%. For a roughened airfoil, the drag coef- 
ficient increases dramatically at high angles of attack because 
of the premature onset of stall. 

Using the results of Brumby, estimates were made of the 
reduction in maximum lift and increase in stall angle for a rain 
roughened airfoil in a landing configuration. 

Using the sand grain roughnesses associated with drop 
impact from Table 8 and thos .ssociated with film waviness from 


Table Ilf the reduction in maximum lift coefficient and in angle 
of attack at stall were calculated. These calculations were 
derived from Brumby's curve for an airfoil with high lift devices 
retracted whose entire upper surface is roughened. 

TABLE 14 

REDUCTION IN MAXIMUM LIFT COEFFICIENT AND 
ANGLE OF ATTACK AT STALL DUE TO ROUGHNESS 


Rain Rate 
(mm/hr) 

Drop Impact 
Cratering 

max 

Film 

Waviness 

(%) 

Drop Impact 
Cratering 

Film 

Waviness 

100 

7% 

11 

1-2* 

1-3* 

200 

13% 

20 

1-3* 

2-4 • 

500 

25% 

25 

2-5* 

2-5* 

1000 

29% 

28 

3-5* 

3-5» 

2000 

34% 

30 

3-6« 

3-5* 


Brumby states however, that leading edge high lift devices even 
in the extended position do not recover degraded lift due to 
large amounts of roughness. 

Table 14 shows a significant reduction of maximum lift and 
stall angle for roughness associated with both drop impact cra- 
tering and film waviness. Though unvalidated assumptions were 
necessary in deriving these results, we believe they are suf- 
ficiently realistic to warrant a most serious consideration of 
the influence of heavy rain on aircraft aerodynamics. Even the 
magnitude of lift penalties and stall angle decreases due to film 
waviness alone-for which our assumptions are best justified— 
could provide serious aerodynamic problems for an aircraft exe- 
cuting a go-around in a heavy rain environment. It is also 
likely that because of the decrease in stall angle, an aircraft 
may actually stall before activation of a stall warning device 
that is based upon clean airfoil aerodynamics. 
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SECTION 9 

AIRCRAFT LANDINGS IN HEAVY RAIN 

We modified an existing digital aircraft landing model 
(Luers and Reeves, 1973) to evaluate the effect of heavy rain 
upon aircraft landings. Each landing begins at 500 feet altitude 
where the aircraft is trimmed by determining the values of angle 
of attack, throttle setting, and elevator deflection, so as to 
result in unaccelerated flight down a 2.7* glide slope. The 
equations of motion were then integrated by a fourth order 
Runge-Kutta scheme. For no rain, the aircraft flies down the 
glide slope at a constant velocity until it reaches the ground. 

An aircraft encountering heavy r?in no longer adheres to the 
glide slope because of momentum and drag penalties due to the 
heavy rain. The impact of heavy rain imparts a reverse thrust to 
an aircraft which decreases its airspeed. The rain roughness 
produces a change in the aircraft drag coefficient. The 
resulting deviation in touchdown point represents the severity 
of heavy rain to landing aircraft. 

Two of the equations of motion are derived by summing the 
forces parallel and perpendicular to V (velocity vector relative 
to the earth) and applying Newton’s Laws of Motion (Luers and 
Reeves, 1973). An additional term in the equation is necessary 
to represent momentum loss due to the heavy rain impact. In the 
parallel direction the term* 

F RAIN ■ costF„ + sintP- 

X X » 

and in the perpendicular direction, the term 

F Rain ■ -cosyF ” sinyF 

2 Z X 

is added to the equation of motion where where y is the flight 
path angle and F^ and F^ are the horizontal and vertical com- 
ponents of the rain momentum. The rain induced drag penalties 


52 











were introduced into the landing simulation by increasing the 
drag coefficient of the aircraft. 





Fixed stick landing simulations were conducted using the 
aerodynamic data for a 747 aircraft and the influence of the rain 
momentum and roughness drag effects were evaluated. No lift 
penalties were introduced. 

Figure 16 shows the trajectories of Boeing 747 landings 
with only momentum effects (Table 2) for five rainfall rates 
(2000inm/hr, lOOOmm/hr, 500mm/hr, 200mm/hr, lOOmm/hr). No drag 
penalty is included in these simulations. Also shown is a fixed 
stick trajectory for a landing in severe wind shear. Without 
rain or wind, the aircraft touches the ground in 10,593 feet. 

With a rain of 2000rom/hr, the shortfall is 2963 feet; for a 
lOOOmm/hr rain, the shortfall is 2051 feet; for a 500mm/hr rain, 
the shortfall is 1211 feet; for a 200mm/hr rain, the shortfall is 
525 feet; and for a lOOmm/hr rain, the shortfall is 255 leet. 

The severe wind shear trajectory resulted from the wind profile 
shown in Figure 17; a strong vertical wind shear of at least 9 
knots/100 feet. The momentum penalty under an incredible rain- 
fall rate does not equal the severe wind shear penalty, neverthe- 
less, the momentum penalty is significant for rainfall rates 
approaching 500mm/hr and above. 


Figure 18 compares the trajectories of Boeing 747 landings 
incorporating only the drag penalties due to heavy rain for five 
rainfall rates with the severe wind shear landing and with a 
landing without rain and wind. The total drag penalties were 
arbitrarily estimated by multiplying the sum of the friction drag 
penalties due to drop impact and waviness by a a factor of 3 to 
account for other sources of drag increase as well as for partial 
roughening on the lower airfoil and fuselage sl rfaces. Without 
rain or wind, the aircraft reaches the ground in 10,593 feet. 

With a rain of 2000mm/hr, the shortfall is 3440 feet, for a 
lOOOmm/hr rain the shortfall is 3226 feet, for a 500ram/hr rain 
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Figure 17, Wind Profile with Large Wind Shear 


Drag Penalty Rain Rate (im/hr) 



Pisure 18. 747 Landings with Drag Penalty Only. Rain Encountered at 300 Feet 


the .horttall is 3000 feet, for a 200mm/hr rain the shortfall is 
2704 feat, and for lOOnm/hr the shortfall is 2245 feet. 

for final comparison, both drag and momentum penalties from 
heavy rain are combined as shown in Figure 19. For a 2000mm/hr 
ra n the shortfall is 3908 feet, for a lOOOmm/hr rain it is 3562 
feet, for a 500mm/hr rain it is 3219 feet, for a 200mm/hr rain it 
IS 2832 feet, and for a 100 mm/hr it is 2343 feet. Under the 
assumed conditions a rainfall rate of approximately 400mm/hr 
would produce the same magnitude penalties of the severe wind 
shear shown in Figure 17. Table 15 summarizes these results. 

In most cases, both heavy rain and severe wind shear will 
accompany each other. With both present, the shortfall is 
greater than with either absent. The increased shortfall is, 
however, not the simple addition of separate shortfalls arising 
from wind shear and from heavy rain. 


LANDINGS WITH GENERAL AVIATION AIRCRAFT 

To our knowledge, no general aviation (GA) aircraft oraches 
h ve been attributed to heavy rain. Heavy rain 1, not pr.clud«3, 
however, aince it haa not been poasible to adeguately identify 
and analyte craah.s in which heavy rain was preeent. in 
addition, GA aircraft have been lesa likely to fly in aeverely 
degraded weather condition, due to pilot and in.trument 
llmitatlona. Aa more and more mnall aircraft become capable of 
natrument flying in poor weather, however, the likelihood for GA 
aircraft encounter, with heavy rain will increaae dramatically. 

Thua, an evaluation of the effect of heavy rain on GA landing, i. 
advisable. 

The aerodynamic charact.rlatlca of a light aingle engine 
high-win, airplane (Greer, et al, 1973) were uaed in almulating 
t e trajectory of a GA plane in heavy rain. The aimulation waa 
aimilar to that for a tran.port aircraft, aa weight, momentum. 
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Rain Rate (m/hr) 



Figure 19. 747 Landings with Drag and Momentum Penalty. Rain Encountered at 300 Feet 







TABLE 15 

SUMMARY OP LANDING SHORTFALLS (FEET) FOR A-747 


Initial Height • 500 feet 
Glideelope ■ 10,593 feet 
Landing Weight - 400,000 lbs. 
Rain encountered at 300 feet. 


Momentum Penalty Only 


Rain Rate (mm/hr) Shortfall 


100 

200 

500 

1000 

2000 


255 

525 

1211 

2051 

2963 


Drag Penalty Only 


Rain Rate (mm/hr) Shortfall 

100 2245 

200 2704 

500 3000 

1000 3226 

2000 3440 


Weight Alone Shortfall 

1% Increase 72 


Drag and Momentum Penalties Combined 


Rain Rate (mm/hr) Shortfall 


100 

200 

500 

1000 

2000 


2343 

2832 

3219 

3562 

3908 
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and drag penalties were separately evaluated under a fixed stick 
assumption. 

The simulation results are similar to those for transport 
aircraft. The weight penalty Is estimated to be less than 3% for 
a small plane; this causes a landing shortfall of less than 200 
feet In a GA landing. The greater maneuverability of small 
aircraft makes such a shortfall negligible. Momentum and drag 
penalties are more severe for GA aircraft than for transport 
aircraft. In all cases, Uie smaller aircraft experiences 
somewhat greater landing shortfalls than transports. 

The fixed stick assumption Is not very realistic for a GA 
aircraft since their slower speeds and quicker responses provide 
greater maneuverability. Much greater latitude la available to 
the GA pilot for remedying penalties Imposed by heavy rain than 
to the transport pilot. Nonetheless, the penalties Imposed by 
heavy rain may create Irrecoverably Inhospitable flylhg con- 
ditions. More attention to realistic scenarios of GA aircraft In 
heavy rain Is warranted. 


SECTION 10 
CONCLUSIONS 

We have investigated some of the aerodynaaic effects of 
rain on an aircraft. The weight penal ty» loss of nomentua due 
to rain impacts » and the increased drag and decreased lift due to 
droplet impact and waterfilm waviness have been considered. The 
weight penalty was found to be insignificant under normal landing 
conditions. The momentum penalty becomes significant for rain- 
fall rates approaching 500mm/hr. Drag and lift penalties could 
be very significant for rainfall rates exceeding lOOmm/hr. 

Landing simulations indicate that the drag and momentum penalties 
alone r associated with a 400 mm/hr rainfall rater may be equiva- 
lent to that caused by a severe windshear of 9 )uiots/100 ft. In 
addition r lift penalties which were not included in the landing 
simulationr would cause even more performance deterioration. 

It was necessary to make a number of approximations during 
the Investigation. A conservative or realistic approach was 
used in contrast to a worst case analysis. For example r in esti- 
mating the momentum penalty of rainr both the deployment of flaps 
and the vertical velocity between plane and air were ignored. 

Both should increase aerodynamic penalties. We believe the esti- 
mates of aerodynamic penalties due to heavy rain# although cruder 
are sufficiently realistic to warrant concerned attention. 

The significance of rain induced momentum loss is con- 
sistent with an antecedent investigation. Evidence for further 
aerodynamic degradation as a result of a rain roughened airfoil 
has been presented here. The combination of the tvo penalties 
have been shown comparable to strong windshear in degrading an 
aircraft landing trajectory. Normally r heavy rain will be accom- 
panied by adverse winds. Separately neither may induce overly 
serious degradation while together both may be lethal. 

A rain-roughened airfoil may cause a reduction of the maxi- 
mum lift coefficient in addition to premature stall. With large 
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amounts of distributed roughness, such as may be caused by heavy 
rain, the deployment of leading edge high lift devices cannot 
reduce this effect. Finally, premature stall may occur insid- 
iously before stall warning devices activate. 

It appears likely that very heavy rain may have been a 
contributing factor in several aircraft accidents. Some of the 
loss in indicated air speed (IAS) appears to have been caused 
by rain. Thus, the derived windshears of accident reconstruc- 
tions may be too large because rain was Ignored. Future research 
will concentrate on refining the effect of rain on drag and 
assessing its effect on lift. In addition, the role of rain in 
several aircraft accidents will be addressed. 
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